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We investigate anomalous damping of the monopole mode of a non-degenerate 3D Bose gas under 
isotropic harmonic confinement as recently reported by the JILA TOP trap experiment [D. S. Lob- 
ser, A. E. S. Barentine, E. A. Cornell, and H. J. Lewandowski (in preparation)]. Given a realistic 
confining potential, we develop a model for studying collective modes that includes the effects of 
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is caused by the joint mechanisms of dephasing and collisional relaxation. Furthermore, the model 
is complimented by Monte Carlo simulations which are in fair agreement with data from the JILA 
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I. INTRODUCTION 

In the late 19th century. Maxwell and Boltzmann un¬ 
covered a path to connect the Newtonian mechanics 
of molecular dynamics to the hydrodynamic equations 
of Euler and Navier-Stokes. These kinetic theories, of 
which the fluid theories were limiting cases, established 
a platform to formulate and analyze the Maxwell dis¬ 
tribution and Demon ma, Boltzmann’s equation and 
the H-theorem, and the assumption of molecular chaos 
[4] (Stosszahlansatz). These fundamental ideas, in par¬ 
ticular the H-theorem, famously stirred controversy in 
the scientific community, which in the 1890s was cen¬ 
tered around the “reversal” paradox and the “recurrence” 
paradox based on Poincare’s theorem [5]. Lesser known, 
but equally as curious was the discovery by Boltzmann of 
a class of exact solutions to his kinetic theory of which the 
Maxwell distribution is a special case [H]. Such solutions 
include the undamped, nonequilibrium oscillations of a 
classical gas under 3D isotropic harmonic confinement: 
the so-called breathing or monopole mode. Whereas 
some of the unexplored aspects in the justification of the 
Boltzmann equation began to be filled in the early 20th 
century, experimental demonstration of the undamped 
monopole mode oscillation has yet to be fully investi¬ 
gated. 

Advances in the trapping of ultracold gases have al¬ 
lowed for the study of collective modes under harmonic 
confinement. Indeed, many experiments have probed the 
transition between the collisionless and hydrodynamic 
regimes [MU in anisotropic traps. Such experiments 
measure the eigenfrequencies and damping rates of the 
various multipole modes and are generally limited by 


three-body losses approaching the hydrodynamic limit. 
Engineering an isotropic 3D trap of sufficiently perfect 
symmetry has remained a technological hurdle in inves¬ 
tigating the undamped oscillation of the monopole mode. 

A recent experiment from the group of Eric Cornell at 
JILA [1] involving a nondegenerate cloud of ®^Rb well 
above the transition temperature utilized a TOP trap 
with anisotropies as low as 0.02%. The monopole mode 
was made to oscillate with minimal damping over many 
trap periods by dithering the trap frequency. However, 
the damping observed was anomalously large given the 
level of anisotropy present in the system. 

This paper explores the inherent limitations on the 
observation of an undamped, nonequilibrium monopole 
mode given a realistic trapping scenario, using the 
JILA experiment as illustration. The effects of slight 
anisotropies in the trap are first quantified and shown 
to lead to coupling between the monopole mode and the 
quadrupole modes. The quadrupole modes suffer from 
collisional damping in the transition between hydrody¬ 
namic and collisionless limits. The importance of cubic 
and quartic anharmonic corrections to the trapping po¬ 
tential are then considered. Such corrections give rise to a 
temperature dependent eigenfrequency and damping rate 
for the collective modes. They also lead to dephasing ef¬ 
fects that create the appearance of actual decoherence 
in the system. In general, such effects generate addi¬ 
tional features observed in the damping in the transition 
between collisionless and hydrodynamic limits. This in 
turn provides an explanation for the anomalous damping 
seen in the JILA experiment. 
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II. BACKGROUND 
A. Boltzmann Equation 

The canonical Boltzmann equation [551 US] describes 
the phase space evolution of a cloud of thermal atoms 
obeying classical statistics: 

+ = /,„;,[/], (1) 


Icou[f] = 
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at a rate determined by the total scattering cross section 
cr. The velocities of the incoming particles are vi and 
V 2 , and the outgoing particles leave with velocities 
and V 2 . For binary elastic collisions, both the center 
of mass (COM) momentum and the magnitude of the 
relative momentum are conserved, and the solid angle 
integral is over the full range of possible final directions 
for the relative momentum. 

As a direct consequence of Boltzmann’s H-theorem, 
the necessary condition for the vanishing of the collision 
integral is that log(/) be composed of quantities that are 
elements in the set of collisional invariants {xi,X 2 -: ■■■}■ 
This defines a class of solutions of the form 

log(/(r,v,t)) = xi+X 2 + •■• (3) 

For a binary elastic collision, this set is spanned by the 
kinetic energy, momentum, and any velocity independent 
function, giving 

log(/(r,v,t)) =a + b-v + cv^, (4) 

where a, b,c are possibly functions of position and/or 
time. The familiar Maxwell-Boltzmann distribution of a 
stationary dilute gas follows from this expression for the 
choice |b| = 0, c = —l/k^T and a = logn . 

Although the parameters that lead to the Maxwell- 
Boltzmann distribution are the most well-known; it 
should be emphasized that Boltzmann noted a curious 
class of parameters that corresponded to collective mo¬ 
tion in a harmonic trap |S]. Boltzmann considered so¬ 
lutions where c had periodic time-dependence describ¬ 
ing temperature oscillations, as well as a position and 
time-dependent b. This class of solutions corresponds 
to the so-called breathing or monopole mode and de¬ 
scribes undamped collective motion of a thermal cloud 
in an isotropic 3D trap, regardless of the oscillation am¬ 
plitude and interaction strength. In a recent paper |5^ . 
the class of undamped nonequilibrium solutions permit¬ 
ted by the H-Theorem was extended to time-dependent 
confinement. For the remainder of the paper only the 
case of harmonic confinement is discussed. 


where /(r,Vi,t) is the single particle phase space dis¬ 
tribution that depends on the generalized coordinate r 
and velocity Vi of all the particles, H is the single par¬ 
ticle Hamiltonian, and {} denotes the classical Poisson 
bracket. The left-hand side describes the single-particle 
evolution of the thermal cloud and the right-hand side 
is the collision integral which drives the system toward a 
state of statistical equilibrium 


|[/(vi')/(v2') -/(Vl)/(V2)] ( 2 ) 


1. Collective modes in a harmonic trap. 


The collective modes of a 3D isotropic trap with trap¬ 
ping frequency Wq reflect the underlying spherical sym¬ 
metry and the quadratic form of the trapping potential. 
Any mode of the trap can be expressed in terms of irre¬ 
ducible spherical tensors, in this case they are multipole 
modes denoted by the labels {l,m). The I = 0 mode is 
the monopole mode of the system, and is characterized 
by an undamped oscillation at 2a;o. The I = 1 dipole 
mode is a COM oscillation independent of the character 
of the interaction and therefore undamped. The I = 2 
quadrupole mode oscillates at 2ujo in the collisionless 
regime and ^/2ujo in the hydrodynamic regime. It is colli- 
sionally damped in the hydrodynamic crossover between 
the two extremes. For small amplitude oscillations about 
equilibrium, one finds a general form for the oscillation 
frequency uj = uJr + iT of the collective mode : 


w 
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^"cL + 


^"hD ~ ^CL 

1 -I- iiOTcoll 


(5) 


where the subscripts denote the mode frequency in the 
hydrodynamic (HD) and collisionless (CL) limits, and 
fcoii is related to TcoU, the thermal relaxation time, by 
Tcoii = {^cl/'^hd)'''coII- The imaginary part of uj cor¬ 
responds to the damping and the real part to the mode 
frequency. This general form is common to many kinds 
of relaxation phenomenon |26j . and, assuming F <C w,. 
can be written in an asymptotic form near the collisional 
and hydrodynamic regimes: 


^ HD 

Fcl 


-A^{ujI^-ujU), 
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( 6 ) 


In the transition regime, again provided that F <C Wr, 
the damping takes on the simple form [26]: 


F 


(7) 
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For traps where the spherical symmetry is broken, the 
monopole mode couples to the quadrupole modes. It 
is therefore collisionally damped at the hydrodynamic 
crossover occurring between the hydrodynamic and col¬ 
lisionless regimes. In Sec. Ill, the impact and strength of 
this coupling on the monopole mode is discussed. 

B. JILA TOP trap experiment ®^Rb 

Of the many schemes for magnetically trapping ultra¬ 
cold neutral atoms, the TOP trap m allows flexibility 
in tailoring the level of anisotropy present in the system. 
TOP traps utilize a fast rotating bias field oriented per¬ 
pendicular to a quadrupole magnetic field. This avoids 
Majorana losses that would otherwise occur due to the 
zero point of the field being inside the cloud [28]. The 
quadrupole and fast rotating bias fields produce a mag¬ 
netic field that has Cartesian components 

/ -I- Bo cos(nf) \ 

B = -fy + Bq sin(flt) (8) 

V -B.z j 


be included in calculations. The potential used in calcu¬ 
lations involving the anharmonic corrections has the form 

u{r, z)ah =f ^^2 + a;2r2) + 

for general coefficients of the cubic and quartic terms (see 
Ref. [29]). 

The monopole mode is driven experimentally by apply¬ 
ing a sinusoidal variation in the trapping frequency over 
four periods of monopole oscillation. The strength of the 
TOP field is changed along with a vertical bias which 
fixes the minimum of the potential (see Ref. [M])- The 
net result is a roughly 25% increase in the mean cloud 
size. 


III. THEORY 

In this Section, two approaches to modeling small am¬ 
plitude collective excitations of a nondegenerate Bose gas 
about its equilibrium state are discussed. 


where Bq is the field strength of the bias field, is the 
strength of the quadrupole field. This bias field rotates at 
angular frequency O, which is typically much larger than 
the Larmor frequency. Near the minimum of the trap the 
cooled atoms experience a time averaged potential 

Uropir, z) = y.Bo + + ?>z^) (9) 

16 Ro 

where /i is the magnetic moment and r = y/ 

In the high-held limit {^lB^ ^ mg) the trap is 
anisotropic with aspect ratio oj^luir = As the 

quadrupole held strength is decreased, the effects of grav¬ 
ity become important as the atoms in the trap sag away 
from the high-held limit trap minimum. The extent 
of the sag is characterized by the dimensionless quan¬ 
tity A = mg/fiBz and shifted trap minimum zq = 
— (BqA) / [BzV^ — A^). The potential expanded about 
the sag position Zq is 

2 ) (1 + A'^)PV1 - 

+ (10) 

The ratio of the trap frequencies is 


LUz 

LOr 
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( 11 ) 


A. Moment method 


By taking moments of the Boltzmann equation Eq. 01 
a set of coupled equations of motion can be derived that 
describe the dynamical evolution of various collective 
modes [22] . For an arbitrary quantity y(r,v), the equa¬ 
tion of motion is 

^-(v-Vr)-/--Vvx\ = (x4o«), (13) 

at \m / 

where the average quantity is given by 

ix) = ^ J rf^?-rf^v/(r,v,t)x(r,v), (14) 


and the collisional average of this quantity is given by 

ixlcou) = ^ J d^rd\iAxIcou[f], (15) 

where Ay = Xi + X 2 - Xi' - X 2 ' with Xi = x(i’. Vi). 
From Eq. (13), it is straightforward to form a closed set 


of coupled equations for the quadratic moments, recalling 
that moments of the form in Eq. Q have no collisional 
contribution. For an isotropic trap, taking the moment 
(r^) yields a closed set of equations describing the evolu¬ 
tion of the monopole mode: 


with A = 7/9 giving an isotropic trapping potential. For 
further details of the modified TOP trap used in the JILA 
experiment, we refer the reader to Ref. |29j . 

Utilizing gravity to symmetrize the trap results in trap¬ 
ping frequencies on the order of 10 Hz. In such a loose 
trap anharmonic corrections become important and must 


d(r^) 

dt 


2(vr)=0, 


d(v • r) 

dt 


^)=0, 


dt 


+ 2ujI{v • r ) = 0 . 


(16) 
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Utilizing a trial solution of the form /(r,v,t) = 


e“‘/(r, v), allows Eq. (16) to be solved algebraically. 


The solution confirms Boltzmann’s general result as pre¬ 
viously discussed in Sec. II for an undamped monopole 
mode at oscillation frequency 2ujo. 

One can also take (Z = 2, m = 0) quadrupole moments 
of Eq. (13) to derive a closed set of six coupled equations 


for an isotropic trap. However, for the quadrupole mode 
the moment, (3u^ — v^) cannot be expressed in the form 
of Eq. ([^, which leads to a nontrivial collisional contribu¬ 
tion {{3vl — v‘^)Icoii) in its equation of motion. For small 
amplitude oscillations about equilibrium, collisional con¬ 
tributions of this form can be approximately rewritten in 
the relaxation time approximation [22l [25l [30] as 


(17) 


iX^COll) - 7 

^coll 

where the relaxation time TcoU is given by 

5 

'^coll — 


^Icoll ' 


(18) 


and 7con = n{0)vthO'/2 is the average collision rate with 
Vth = \/SkbT /ttto. As with the monopole mode, the 
set of coupled equations for the quadrupole mode can 
also be solved algebraically, giving frequencies of 2a;o and 
\/2wo in the collisionless and hydrodynamic regimes re¬ 
spectively. However, in the hydrodynamic crossover the 
quadrupole mode damps as in Eq. ([^ due to the nonva¬ 
nishing collisional contribution. 


B. Scaling ansatz 


The scaling ansatz method has been used for Bose 
gases above and below the transition temperature [3T]- 
133]. In the nondegenerate regime, the ansatz method 
has been used to evaluate the effects of the mean field 
interaction [331 [33] and to estimate the effects of anhar- 
monic corrections to the trapping potential on the collec¬ 
tive mode frequencies and damping rates |34| . Here this 
method is outlined. 

When the thermal cloud is in statistical equilibrium 
under harmonic confinement, the Boltzmann equation for 
the equilibrium distribution /o(r,v) is just the Poisson 
bracket: 


9/o 2 9fo \ 

=0, 


dv. 


(19) 


where the subscript denotes components of a Cartesian 
vector in 3D and the possibly differing frequency com¬ 
ponents in each direction. Moments of the equilibrium 
distribution are the same as those obtained in the previ¬ 
ous subsection when the temporal derivatives are set to 
zero. 

To describe time-dependent collective oscillations, a 
scaling ansatz on the form of the distribution function 


/(r,v,t) can be made: 

/(uv,t) = r/o(R(r,t), V(r,v,t)). (20) 

This ansatz utilizes the symmetry of the problem and 
the nature of the collective modes in the forms of the 
renormalized position, R(r, t), and velocity, V(r,v,t), 
with an additional factor T to enforce the normalization. 

For the dipole mode, the form of the ansatz is Ri = r^ — 
r]i(t) and Vi = Vi—rii{t) with T equal to unity. The details 
of the motion are contained in the time-dependent vector 
of free parameters, ff. Such a scaling clearly describes a 
translation of the COM of the equilibrium cloud. 

For the monopole and quadrupole mode, the ansatz 
mimics the form of Boltzmann’s solution as discussed in 
Sec. H: 





( 21 ) 


where there are two time-dependent vectors of free pa¬ 
rameters 9 and b. Such a scaling describes the stretching 
and compression of the equilibrium cloud in phase space 
along with a space-dependent and time-dependent trans¬ 
lation in the local velocity, while maintaining a stationary 
COM. 

Substituting Eq. (20) into Eq. ([^ and following [311 
133] . a Newton-like set of equations of motion for the free 
parameters can be derived 


k + 2 ^ 0 , =- X [0^-0]^ ( 22 ) 

'^coll 

where the quantity 9 = J2i is the average tempera¬ 
ture. To obtain information about the collective oscilla¬ 
tions, the free parameters are linearized about the equi¬ 
librium position for small amplitude oscillation, « 1 
and 0i « 1, and assumed to have time-dependence of 
the form For a harmonic trap, the results obtained 
agree with the moment method. 


C. Modeling trap imperfections 


Slight anisotropies in the harmonic confinement and 
anharmonic corrections to the trapping potential Eq. (12) 


present possible sources of damping. Modeling these im¬ 
perfections is discussed in this subsection using trapping 
data from the JILA experiment. 
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FIG. 1. Monopole damping rates versus collision rate for 
a set of frequencies from the JILA experiment with various 
anisotropy parameters and central frequencies fo- The legend 
gives the level of anisotropy, where A, e, and (Uq correspond 
to rewriting the harmonic confining potential as in Eq. (23l. 


1. Damping due to trap anisotropies 

Anisotropies in the confining potential lead to coupling 
between the monopole and quadrupole modes, which are 
damped in the hydrodynamic crossover regime. The mo¬ 
ment method can account for this coupling in a straight¬ 
forward way given a potential of the form 

U{x, y,z) = ^ ((1 + + (1 - e)y^ + X^z^), (23) 

where A = oJzIujq, (1 -I- e)wo = oJx, and (1 — e)uJo = tOy. 

Starting with (r^), the moment method yields a set of 
nine coupled equations for the monopole and quadrupole 
modes (see Ref. [35]). In Fig. 1, frequency measurements 
from the JILA experiment with various anisotropies were 
used to estimate the damping rate for small oscillations. 
For nonzero e and A, the monopole couples to multi¬ 
ple quadrupole modes and the shape of the curve de¬ 
parts from that expected for the typical transition region 
shape when e = 0. The monopole damping measured 
in the JILA experiment is on the order of 0.1 s“^ or 
larger. Therefore, the trap can roughly be treated as 
isotropic with A = 1 and e = 0 for the lowest four curves 
in Fig. 1. Thus, the observed monopole damping must 
then be due to some other mechanism that is present 
even when the harmonic piece of the trapping potential 
is virtually isotropic. 


2. Damping due to anharmonic corrections 

Computing the monopole damping using the moment 
method for nonzero anharmonic coefficients generates an 
infinite set of coupled equations for successively higher 
order moments and so is not a viable method. However, 


the scaling ansatz method, which takes as an input mo¬ 
ments of /o, can provide an estimate for the effects of 
the anharmonic coefficients which enter through higher 
order moments of /q. The usefulness of such an approach 
hinges on being able to neglect deformations of the equi¬ 
librium distribution due to the anharmonic corrections. 


Using the scaling ansatz, Eq. (22) can be rewritten in 


a more general form for an arbitrary trapping force, Fi: 

h - ^{F^ibJrj)ri)o = 0 , 

0 , + 2 ^ 0 , =- —[ 0 ,- 0 \, ( 24 ) 


where the 0 subscript indicates a moment as in Eq. (14) 
over the equilibrium distribution /o, and bjVj is a 
short hand notation for the set {bxrxibyry^bzVz}- Erom 
Eq. (24), the procedure for including even order correc¬ 


tions to the trapping potential is straightforward. How¬ 
ever, to include odd order corrections, the equilibrium 
distribution, /o, must be deformed. A simple deforma¬ 
tion is 


fo 


ah 


(r,v)= 1- 


mj3 


(x^ +y‘^)z- 


ma 

3fchT' 


/o(i-,v) 


which is a perturbative expansion to first order in the 
odd order anharmonic corrections from Eq. (12). The 


deformation does not change the overall normalization 
nor the value of the even-order moments, and the added 
terms are collisional invariants of the form Eq. (|^. 

In Sec. VI, the validity of this deformation compared 
to numerical and experimental results for the damping 
of the monopole mode is addressed. However, in the 


remainder of this Section, the general solution of Eq. (24) 


is discussed for small amplitude oscillations, following the 
derivation in Ref. |34] and noting that the derivation is 
the same for /g^. 

Eq. (24) is solved in both the collisionless and hy¬ 
drodynamic regimes. Making the assumption that the 
monopole damping is much less than the collective mode 
frequency (Fm ^ ujm) leads to a generally applicable for¬ 
mula Eq. 0 for the damping. In the JILA experiment, 
^m/ojm is on the order of 10“^ to 10“^. 

In the collisionless regime, TcoU —t oo, which gives a 
power law relation between the two scalin g pa rameters: 


A. Using this power law relation, Eq. (|24|) collapses 


into three differential equations for each component of 
the parameter bi. Eor small amplitude oscillations bi « 
1, and the substitution rji = bi — 1 linearizes the three 
differential equations for each component of yf. 


.. , 3 (C/iri)o 

"I-r“2^- 


m 


(uAo 


1 


m 


3 


{Utjnrj)o 

(r!)o 


Vj = 0, (26) 


where Ut = dU/dri and Uij = d'^U/dridvj. 
be treated with matrix methods as in Ref. j 
the collective mode frequencies. 


Eq. (26) can 


to extract 
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In the hydrodynamic regime, TcoU —>■ 0, and thus the 
system is in local equilibrium everywhere. Furthermore, 
the components of the 6iS must be equal to to their av¬ 
erage, 6i = 0. From the normalization of the ansatz 
then follows a power law relation between the parame¬ 
ters: 9i — 1/ rij ■ This implies that Eq. (24) collapses 
to a set of three differential equations for eacncomponent 
of the parameter hi that couple when the power law re¬ 
lation for 9i is substituted. Making the small amplitude 
assumption leads to the following equation: 


1 {U%ri)o 
m (r2)o 


{Uzjrirj)o 



(27) 

As in the collisionless case, matrix methods can be 
used to extract the collective mode frequencies for the 
monopole mode and the (1=2,m=0) and (1=2,m=2) 
quadrupole modes. 


The calculated mode frequencies have temperature de¬ 
pendence through higher-order moments coming from 
the anharmonic corrections. For a given temperature 
there is a shift in the frequency and damping from that 
anticipated by the harmonic limit. In the collisionless 
limit, this shift in frequency gives rise to dephasing- 
induced damping and creates the appearance of actual 
relaxation in the system over short time scales. As dis¬ 
cussed in [31], it is expected that in the collisionless 
regime the dephasing induced damping should scale pro¬ 
portionately to the shift in the mode frequency. Aw « 
AF. In general, it is less clear how Aw correlates with 
the shift in the damping rate. Therefore a transition 
regime exists where dephasing effects in the cloud are 
destroyed through collisions, and collisional damping be¬ 
comes the dominant effect. Such an effect can account 
for anomalous damping of the monopole mode in regimes 
where ujmTcoU S> 1, whereas collisional damping dom¬ 
inates when lumTcoU ~ 1- We refer to this regime as 
the dephasing crossover, and in Sec. V this crossover is 
studied in detail with the aid of numerical simulation. 


IV. NUMERICAL METHODS 


In this section the numerical algorithm used to simu¬ 
late the dynamical evolution of the thermal cloud is de¬ 
scribed. Numerical simulation of the system isolates the 
roles of various damping mechanisms and permits quan¬ 
tification of the dephasing effects in the cloud. The ther¬ 
mal cloud algorithm from Ref. |36j is adopted, working in 
the quantum collisional regime m where two-body col¬ 
lisions are s-wave. However, the many-particle statistics 
are still classical and well-described by the Boltzmann 
equation. The simulation consists of a swarm of tracer 
particles that act as a coarse-grained distribution func¬ 


tion [55] : 


/(r, P, t) « ^ '5(1' - ri)<5(p - Pi), (28) 

i 

where Ntt/Ntp is a weighting factor. The sum is over 
the entire set of tracer particles, where each is uniquely 
described by their position and momentum (ri,pi). The 
tracer particles first undergo collisionless evolution via a 
second order sympletic integrator [39l |40| . Following free 
evolution, the tracer particles are binned in space and 
tested for collisions. 


A. Collisions 


The collision algorithm, which is a procedure for Monte 
Carlo sampling the collision integral m, follows that de¬ 
scribed in Ref. [SS]. After binning the particles in space, 
random pairs in each bin are selected and their collision 
probability is calculated. If a collision is successful, the 
particle velocities are updated according to the differen¬ 
tial cross-section for an s-wave collision. 

For each set of simulation parameters we check that the 
equilibrium collision rate averaged over the entire cloud 
matches with the analytic result: 


Teg — 


auj^m 
2Tr^ kbT’ 


(29) 


to < 2%. The velocity independent cross-section used is 
that for bosons a = Stto^. 


B. Drive mechanism 

To mimic the experimental monopole drive scheme, 
the trap frequency, wq, is modulated over four periods 
of monopole oscillation 

where A is a unitless measure of the strength of the 
drive. During the drive, the anharmonic corrections are 
neglected, and are switched on when the next oscillation 
minimum occurs after the drive is turned off. Fig. 2 illus¬ 
trates the increase of the monopole amplitude during the 
drive and the effect of turning on the anharmonic shifts 
on the total cloud energy. For A = 0.15 and a starting 
temperature of 152nK, the percent increase in the to¬ 
tal energy when the anharmonic corrections are switched 
on is on the order of 0.001%. During the drive there is 
also a shift in the average temperature and full width at 
half maximum (FWHM) of the oscillation due to an in¬ 
crease in the energy of the cloud from the work done on 
it by pumping of the trap. It would be prohibitive to in¬ 
clude the modulation of the bias, quadrupole, and all of 
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FIG. 2. Driving the monopole moment for A = 0.15 and ini¬ 
tial temperature of 152nK over four periods of oscillation, (a) 
FWHM of the monopole moment versus holding time. The 
grey area indicates the time in between the end of the drive 
and switching on the anharmonic terms in the trapping po¬ 
tential. (b) The % change in the total energy of the cloud 
between a trapping potential that is purely harmonic and the 
same potential with anharmonic terms switched on. The ver¬ 
tical line indicates the point at which the anharmonic terms 
are added to the trapping potential when the % increase of 
the total energy is on the order of 0.001%. 


the shimming fields used experimentally to recreate the 
drive. The simulation driving scheme effectively captures 
the experimental drive. Results are presented then over 
a range of drive strengths which produce an increase in 
the mean cloud side at the end of the drive phase which 
is on the order with that experimentally measured. 


V. RESULTS 


We are now in a position to compare and contrast the 
results of the JILA experiment against the theoretical 
model and numerical simulation. Here, only data with 
the monopole drive is considered. Moreover, our analysis 
is focused on experimental data that is not dominated by 
effects due to trap anisotropy, as this is well-understood 
[13 ISO]. Characterizing the sensitivity of the monopole 
mode to anharmonic corrections around the dephasing 
crossover is the main result of this paper. This sensitiv¬ 
ity stands as a general issue for undamped, nonequilib¬ 
rium collective modes, such as the monopole oscillation. 
The quadrupole modes must also damp around the de¬ 
phasing crossover; however, these modes are not explored 
further as they do not fall into the category of undamped, 
nonequilibrium modes. 

In the JILA experiment, the monopole data was taken 
for a range of atom numbers and temperatures. The 
small cloud data {N « 10^ atoms) which lies in the col¬ 


lisionless regime is analyzed first. Finally the large cloud 
data {N « 10^ — 10® atoms) which lies between the col¬ 
lisionless and hydrodynamic regimes is analyzed. 


A. Collisionless regime 


In the collisionless regime, the many-body dynamics 
of the thermal cloud are dictated by the single-particle 
trajectories. Here, the system behaves in analogy with 
a simple pendulum. For small amplitude oscillations the 
pendulum executes simple-harmonic motion. However, 
as the amplitude grows, the small-angle approximation 
breaks down and anharmonic corrections become impor¬ 
tant. The pendulum traces out a repeating trajectory 
but with an energy-dependent period. 

The procedure for replacing anharmonic corrections by 
an energy-dependent harmonic potential in ID is well- 
known [2bj . The renormalized trapping potential for 
small amplitude motion along the z-axis is: 


U{0,0,z) 

CJz 




mQl{E) 2 
2 ^ ’ 

^z(l + ^zE), 

5a^ 3k 


(31) 


In Fig. 3(a), atoms in the anharmonic trap have been 
binned according to their energy and the collective mo¬ 
tion in each bin averaged to obtain the monopole period 
in each bin. For the first several bins the period of col¬ 
lective oscillation increases linearly with the bin energy. 
For higher energy bins the scaling is still monotonic; how¬ 
ever, the deviation increases as the difference in the linear 
correction between the different axes becomes more ap¬ 
parent. At sufficiently high binning energy, the linear 
correction Eq. (31) breaks down as higher-order terms 
become important, and the atom number in each bin 
decreases, spoiling the appearance of an undamped col¬ 
lective mode. 

In the absence of collisions, the population of atoms 
in each bin is static over the entire simulation. The sta¬ 
tistical weight of each frequency component, as shown 
in the normalized histogram of Fig. 3(b), is then also 
static. The determining characteristics of the frequency 
distribution are its width and shift of the peak from 
= 2tLio. The theoretical prediction of Alum from the 
scaling ansatz using is for the shift of the peak (grey 
region of Fig. 3(b)) and agrees with the numerical re¬ 
sult. However, the width, denoted Sojm, of the frequency 
distribution determines the relative importance of de¬ 
phasing effects. If the frequency distribution is shifted 
and sharply peaked, dephasing induced damping is min¬ 
imal, mimicking the large amplitude oscillation of a sin¬ 
gle pendulum. Otherwise, if the distribution is broad, as 
in Fig. 3(a), the net sum of many ‘pendulums’ oscillat¬ 
ing at shifted frequencies is dephasing induced damping 
through interference. 
















Bin energy (fc(,r) 



0 0.01 0.02 0.03 

1 - uj/u)li 


FIG. 3. Resolving the energy dependence of the monopole 
moment for a drive strength of = 0.15 (9 ± 1.35Hz) and 
an initial temperature of 152nK. (a) After binning the tracer 
particles into energy bins, each bin oscillates as an indepen¬ 
dent monopole. Here, the ratio of the average period of 
the monopole in each bin and the zero temperature result, 
T ^4 = along with the standard deviation (grey re¬ 

gion), is plotted versus the energy of each bin in units of 
thermal energy, (b) The statistical weight, P{uj), of each 
frequency component in the cloud. The grey region is the 
theoretical prediction for the shift of the monopole frequency 
from L (150nK) to R (ISOnK). 



Drive strength 



We now compare directly to the small cloud data from 
the JILA experiment. Starting with the temperature 
and atom number quoted in the experiment, the drive 
strength is varied and the amplitude of the resulting os¬ 
cillation is fit to a simple exponential decay to compare 
directly to the fitting function used in the experiment. As 
shown in Fig. 4(a), a drive strength of A = 0.15 matches 
the experimental data, suggesting that the anomalous 
damping seen in the experiment for the small cloud data 
is mainly due to dephasing effects. Such effects depend 
on the settled temperature of the cloud as shown in 
Fig. 4(b). 

Although the width of the frequency distribution de¬ 
termines the strength of the dephasing induced damping, 
in the collisionless regime where the distribution is broad, 
one naively expects that the frequency shift Awm and the 
width scale proportionately and can be used as a fitting 
function, looking then for a scaling fit of the form 


FIG. 4. (a) Damping of the monopole over a range of drive 
strengths A = {0.05, 0.1, 0.15, 0.20,0.25} and temperatures in 
the collisionless regime. The grey regions are the damping 
results along with the uncertainty from the JILA experiment 
for (from bottom to top) 65nK, 125nK, and 152nK. (b) The 
dashed lines are the theoretical predictions for Aujm{T) using 
/o ^ (single dash) and using /o (double dash) for the scaling 
ansatz versus the settled temperature of the cloud. The points 
are the damping results from (a). The solid line is a fit to the 
numerical damping results using Eq. (321. 


Fig. 4(b), as the settled temperature increases the de¬ 
formed gaussian predicts that the peak frequency de¬ 
creases {Aujm{T) > 0) as shown by the dashed lines in 
Fig. 4(b). 


B. Crossover regime 


^dephase{T) = Awm(C^), 


(32) 


where C is a free parameter independent of the temper¬ 
ature. For the experimental trap parameters, C = 0-4 
from a fit to the data shown in Fig. 4(b). 

The scaling ansatz result using /o instead of /q ^ pre¬ 
dicts a decrease of the period with increasing bin en¬ 
ergy, which disagrees directly with the numerical result 
and Eq. (31). Additionally, from the dashed lines in 


The collisionless regime picture of a collection of un¬ 
coupled monopole modes oscillating at shifted frequen¬ 
cies begins to break down when the dephasing period 
Tdph = 2Trj6oJM exceeds TcoU- This is typically in a dif¬ 
ferent collisional regime than the hydrodynamic crossover 
where ujqTcoU ~ 1, and an average atom suffers a col¬ 
lision on a timescale faster than the oscillation period. 
Therefore, as the collision rate increases there are two 
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FIG. 5. Damping rate versus the equilibrium collision rate. 
The triangle data points are a result of numerical particle 
simulations. The blue curves are scaling ansatz predictions 
using /o ^ for collisional damping of the monopole mode for 
temperatures 60nK, 80nK, and lOOnK from bottom to top, 
respectively. 


important regimes: the dephasing crossover and the hy¬ 
drodynamic crossover. 

SujmTcoU ~ 1 (Dephasing Crossover), (33) 

^qTcoU ~ 1 (Hydrodynamic Crossover). (34) 


In the remainder of this subsection an individual ex¬ 
perimental run from the large cloud data set at 60.6nK 
with 3.463 x 10® atoms with measured trap frequencies 
ifx = 9.036Hz, fy = 9.034Hz, = 9.034Hz) is analyzed. 

From the criteria Eq. (341 can be obtained crude esti¬ 
mates for the dephasing crossover jcoii ~ 10“^ — 1 
and the hydrodynamic crossover 'ycoU ~ 10^ s“^. This 
data set corresponds to the lowest curve in Fig. 1 for a 
quoted collision rate 'ycoii = 8-88 s“^ which lies between 
the two regimes. The maximum damping rate due to 
anisotropies from Fig. 1 is on the order of 10“® s“^ com¬ 
pared to the experimentally measured rate 0.14 ± 0.02 
s“®. Therefore, the harmonic part of the trapping po¬ 
tential is treated as isotropic. 

To quantify the effect of the anharmonic corrections 
in the crossover regime, numerical simulations with A = 
0.05 over a range of collision rates using the trapping 
frequency /o = 9.035Hz and initial temperature 60nK 
were performed. Fig. 5 contains the results along with 
the prediction for the collisional damping from the scaling 
ansatz theory using /q^. The scaling ansatz prediction 
is plotted for several temperatures, illustrating the effect 
of the increase in temperature of the cloud post-drive on 
the damping. 

The estimates for the hydrodynamic and dephasing 
crossover regimes also agree with the qualitative struc¬ 
ture of Fig. 5. The dephasing crossover is marked by 
a decrease in the dephasing induced damping; whereas 
the hydrodynamic crossover is marked by an increase in 
the collisional damping. The region in between the two 



Bin energy (fc^T) 


FIG. 6. Resolving the energy dependence of the monopole 
moment for A = 0.05 and initial temperature of 60.6nK for 
a range of total atom numbers. Each data point is the peak 
component of the Fourier transform of the normalized collec¬ 
tive oscillation in the bin. From bottom to top the curves 
represent frequency spreads in the collisionless, dephasing 
crossover, and hydrodynamic crossover respectively. 


crossovers is characterized by a local minimum in the 
damping rate. Comparing the 65nK damping rate from 
Fig. 4 with the 60.6nK damping rate from this subsection 
draws experimental support for this result. That the col¬ 
lisionless regime is characterized by higher damping than 
the dephasing crossover is a common feature of all of the 
large cloud data compared for similar temperatures to 
the small data [T]. The intermediate regime between de¬ 
phasing and hydrodynamic crossovers then provides an 
experimental window where damping from trap anhar- 
monicities can be minimized. 

Whereas the physics of dephasing in the collisionless 
regime is analogous to the uncoupled oscillation of a col¬ 
lection of pendulums with different oscillation energies, 
in the dephasing crossover the pendulums begin to cou¬ 
ple, and, as illustrated in Fig. 6, the width of the spread 
in frequency components narrows. The coupling effec¬ 
tively synchronizes the different pendulums, which leads 
to a local minima in the damping rate as seen in Fig. 5. 
As the coupling increases, the system nears the hydro- 
dynamic crossover, and the spread in the frequencies is 
minimal. The damping is then due to the appearance 
of a temperature dependent anisotropic trap, where the 
monopole mode damps naturally through coupling with 
the quadrupole modes whose moments are not collision- 
ally invariant. 


VI. CONCLUSION 

In a 3D isotropic trap the monopole mode executes 
undamped nonequilibrium oscillations as predicted by 
Boltzmann. Such a setup has been recently realized 














10 


experimentally at JILA experiment; however, this ex¬ 
periment measured anomalous damping rates for the 
monopole mode. In this paper, possible sources of damp¬ 
ing given a realistic trapping scenario based on the JILA 
experiment are presented. Slight anisotropies in the trap¬ 
ping frequencies were shown to cause negligible damping. 
In contrast, anharmonic corrections to the trapping po¬ 
tential can account for the damping and persist even in 
the collisionless regime where dephasing effects mimic ac¬ 
tual decoherence of the system. As the collision rate is 
ramped up, the system traverses the dephasing crossover 
regime, which is characterized by a local minimum in 
the damping rate. A reduction in the damping is consis¬ 
tent with the rates observed experimentally and provides 
a region of minimum damping where the effects of trap 
imperfections are reduced. Quantifying the impact of an¬ 


harmonic corrections is critical to characterizing damping 
in the monopole mode in virtually isotropic traps. 
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